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Abstract 

The statistical modelling of spatial extremes has recently made major advances. Much of its focus so far has been 
on the modelling of the magnitudes of extreme events but little attention has been paid on the timing of extremes. To 
address this gap, this paper introduces the notion of extremal concurrence. Suppose that one measures precipitation 
at several synoptic stations over multiple days. We say that extremes are concurrent if the maximum precipitation 
over time at each station is achieved simultaneously, e.g., on a single day. Under general conditions, we show 
that the hnite sample concurrence probability converges to an asymptotic quantity, deemed extremal concurrence 
probability. Using Palm calculus, we establish general expressions for the extremal concurrence probability through 
the max-stable process emerging in the limit of the componentwise maxima of the sample. Explicit forms of the 
extremal concurrence probabilities are obtained for various max-stable models and several estimators are introduced. 
In particular, we prove that the pairwise extremal concurrence probability for max-stable vectors is precisely equal 
to the Kendall’s t. The estimators are evaluated by using simulations and applied to study the concurrence patterns 
of temperature extremes in the United States. The results demonstrate that concurrence probability can provide a 
powerful new perspective and tools for the analysis of the spatial structure and impact of extremes. 
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1 Introduction 

While most of the time extreme value analysis focuses on the magnitude of extreme events, i.e., how large extremes 
events are, little interest has been paid to their genesis. This paper tries to fill in this gap by looking at what we shall call 
concurrency of extremes, e.g., have two locations been impacted by the same extreme event or was it a consequence 
of two different ones? Eor example, one could observe daily rainfall at various weather stations and would like to 
quantify the risk that the rainfall extremes over a spatial domain are due to a single extreme event, i.e., a large storm, 
affecting the entire area. Although potentially rare, such events have great socio-economic consequences and their 
probabilities should be assessed precisely. 

More formally, given a sequence Ai,..., A„ of independent copies of a stochastic process X defined on a compact 
set X c R"', d> 1, we say that extremes are sample concurrent at locations si,..., s^ e X, k > 2, if 

max Xiisj) = X{isj), j = l,...,k, (1) 

for some { e Clearly this means that only the observation X( contributes to the pointwise maxima at 

locations si,.. .,Sk. It occurs with probability 


Pn{Sl,...,Sk) =P 


for some ^ e {1,; 


max Xiisj) = X{(sj), j ^l,...,k 


( 2 ) 
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henceforth referred to as sample concurrence probability. 

Provided that X has continuous margins it is not difficult to see that 

p„(_Sx, . . . , Sk) = «E 

where F is the multivariate cumulative distribution of .. .,X{sk)}. Interestingly, concurrence of extremes event 

is invariant under increasing transformations of the marginals so that the concurrence probability does not depend on 
the marginal distributions of X but only on its dependence structure, i.e., the copula C associated to F. 

One drawback of the sample concurrence probability p„(si,..., s*) is that it varies with the number of observations 
n. Surprisingly, however, we prove in Theorem[T]below that, under mild regularity conditions, this quantity stabilizes 
to a universal large sample limit 


p„{si, ...,Sk) —> p{si,Sk), n OO. 


(3) 


Throughout this paper, we will call the limiting probability p(si, ..., Sk) the extremal concurrence probability. This 
asymptotic quantity is naturally expressed in terms of a max-stable process 77 emerging in the limit of the normalized 
maxima in O, as 77 —> OO. A direct, intuitive, and equivalent definition of the extremal concurrence probability can be 
given in terms of the spectral r epresentation of this max-stable proce ss rj. 

Following ide Haan 1 1984 1. Penrose 1 1992l] and Schlatherl 1 2002 1. let 


77 ( 5 ) = max 7 ,( 5 ) s 6 A, (4) 

/>1 

where ; / > 1) are the points of a Poisson process on (0, 00 ) with intensity measure T, are independent 

copies of a non negative stochastic process with continuous sample paths such that E{T(i)) = 1 for all i 6 A and 
E{sup^g_;^- T(i)) < OO. It is often more convenient to rewrite (HI into 

77 ( 5 ) = max ^(i), s e X, (5) 

(^g<D 

where O = {i^o,; i > 1) with (p, - (jYi is a Poisson point process on Co, the space of non negative continuous functions 
on X. Within this framework, we now say that extremes are concurrent at si,..., sj: e A if 


T](sj) ^ ipfisj), j^l,...,k, (6) 

for some { > 1, and similarly to the definition of the sample concurrence probability the extremal concurrence 
probability is defined by 


p{si, ...,Sk) = P {for some (>1: rj(sj) = (fidsj), 7 = 1,.. .,k}. (7) 



6»(ii,i2) = -l0gP{77(ii) < 1 , 77 ( 52 ) < 1), 51,52 6^- (8) 


Interestingly, the extrem al concurrenc e probability and the pairwise extremal coefficient share connections. For in¬ 
stance Proposition 5.1 in IStoevI 020081] implies 


^{2 - 0(5i, 52 )) < pis\,S 2 ) < 2{2 ■ 




(9) 


and we shall see later that the properties of the extremal concurrence probability are similar to that of the pairwise 
extremal coefficient. 

The structure of the paper is as follows. In Section|2]we make connections between sample concurrence probabil¬ 
ities and their extremal counterparts and derive their properties. Section |3] gives closed forms for various parametric 
max-stable models, and Section|4]introduces various estimators for the sample/extremal concurrence probabilities. The 
proposed estimators are then analyzed in a simulation study in Section |5] and applied to US continental temperature 
extremes in Section|6l 
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2 Concurrence of extremes 


In this section we show that sample concurrence probabilities converge to extremal concurrence ones under rather 
mild domain of attraction conditions. We then provide formulas for the extremal concurrence probability based on the 
spectral representation of the associated max-stable process and establish their basic properties. 


2.1 Sample and extremal concurrence 


Concurrence of extremes can be defined through the more general notion of a hitting scenario, which reflects precisely 
how many different events contribute to the componentwise maximum. Let Xi,... ,X„ be a sequence of independent 
copies of a stochastic process X defined on X and si,.. .,sic 6 .Y be different locations. We suppose that X has 
continuous marginals to ensure that {Xi(iy),.. .,X„{sj)} has no ties almost surely and that the maximum is uniquely 
reached. Let Mn(s) - max,=i . X,(i) be the componentwise maximum and consider the sets C, = {j: M„(sj) = 
Xi{sj)}, / = that account for the location where the f-th component X, dominates the rest. Some of these sets 

may be empty, but from the above discussion, with probability one the non-empty ones are disjoint and form a random 
partition of {1,..., A:). This partition n,, - {C,: C, ^ 0} will be referred to as the sample hitting scenario. 

By analogy with extremal concurrence, one ca n define an extremal hitting scenario associated to a max-stable pro¬ 
cess b y using the underlying Poisson point process | Wang and Stoev , 2011 ; Dombrv et al.L 20131 Dombrv and Evi-Minko 
201311 . More precisely, for 77 as in (|5]l, the extremal hitting scenario n is defined as the random partition of {1,..., A:) 
such that two indices 71,72 6 {1,..., A:) are in the same component of n if and only if 


argmax = argmax ipi(sj^). 

i>l i>l 

Whatever type of concurrence is considered, i.e., sample concurrence ([TJ or extremal concurrence (|6l), extremes are 
said concurrent if and only if 7r„ = {1,..., A:) or tt = {1,..., A:). The next theorem shows the convergence of the sample 
hitting scenario to the extremal one. 


Theorem 1. Assume that [/ii{X(ii)),..., hk{X{sk)W belongs to the maximum domain of attraction of{T]{si ),..., T](sk)}, 
for some strictly increasing deterministic functions hi, i = 1,... ,A:. Then, the sample hitting scenario 7 t„ converges 
weakly as n ^ 00 to the extremal hitting scenario n associated to the max-stable process rf. 

Proof Let Xi,X 2 ,... be a sequence of independent copies of X observed at some locations i = (ii,..., ii:) 6 X^, 
k > 2. For brevity we will use componentwise algebra (sum, maximum, etc.) and write X(s) = {X(si),... ,X(sj:)). 
Since the hitting scenario is invariant to strictly increasing deterministic transformations of the marginals, we can 
assume without loss of generality that hi(x) - x, i = 1,... ,k. By assumptions the componentwise maxima M„{s) = 
{max/=i „ Xfsi ),..., max/=i „ converge in distribution 


Mnis) - b„{s) 
a„(s) 


his). 


n —> 00 , 


where a„(i) > 0 and b„(s) 6 R. It is well known that the above convergence is equivalent to the convergence of the 
point process 


o„ = 


Xfs) - bnis) . 


a„{s) 


: i - 1,... ,n> (Z 


to a Poisson point process O on [0, oo)^\{0) where the convergence is meant in the space of point measures ^^([0, oo]*\ 
{0)) equipped with the metric of vague convergence llResnickLll987 l. 

Consider the mapping 

p([0,oof\{0})^lf^k 
T' n('P), 


n: 


where is the set of all possible partitions of {1,..., A:) and n('P) is the hitting scenario associated to the collection of 
functions T* = { 1 / 7 ,: / 6 /), / c N. Clearly, 7r,j = n(<l)„) and n - n(<l)). Since the map 11 is well defined and continuous 
at each point 'P 6 .y£p{[<d, 00 ]*^ \ (0)) for which the maxima are uniquely defined, see Apnendix lA.ll we can apply the 
continuous mapping theorem to show that the weak convergence <1>„ —> O entails the weak convergence 7r„ —> tt. □ 
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Remark. A direct consequence of Theorem[T]is that, provided the stochastic process X is in the domain of attraction 
of some max-stable process 77, the sample concurrence probability converges to its extremal counterpart, i.e., 

Pn{s\,...,Sk) = P[7r„ = {1,...,^)] —» P[7r = {1,...,^)] = p{si,...,Sk), n ^ 00, 

and proves Q. 




2.2 General formulas and properties of extremal concurrence probabilities 

The following theorem gives an expression for the extremal concurrence probability p(si ,..., Sk). 

Theorem 2. 'We have 

i\ ( Y(s:) 

p{si,. .., Sk) = Ey Ey i max 

[[ 1^7=1,...,7: Y(Sj) 

where Y and Y are independent copies of the stochastic process appearing in (|4]l. 

Proof. Let A denotes the intensity mesaure of the Co-valued Poisson point process (t> in Q given by 


A(A) = r P(^T e A)r^d^, 

Jo 

for all Borel set A c Co- We have 

p(si,...,sk) = P{39? 6 cp: 95(^1) = 77(ii), . ..,(p(sk) = T](sk)} 
¥[f(sj)>fi(sj)J^l,...,k}A(df) 

r 


fl 

*JCo 


r 


exp 


= Ey 


-Ef I max 


Y(sj) 
f=C'..k ^Y(Sj) 




X max - 


( 10 ) 


( 11 ) 


where ij and Y are independent copies of 77 and Y respectively. Note that the second equality uses Slyvniak’s formula, 
the fourth one the cumulative distribution of the max-stable process 77 and the last one the expectation of an inverse 
exponential random variable. □ 


Remark. By the seminal paper of de HaanI 1 1984 ] (see also Stoev and Taaau 1 2005 1 ; Kabluchl^ 1 2009l] ). any continu¬ 
ous in probability max-stable process can be represented as 


{ri(s)\ i 6 A) = (max^,/,(M/): 
1 i>l 


6 A 


( 12 ) 


where I/,: s 6 A) is a collection of non-negative integrable functions on the space (U,'L(, v). Here {(^,, 17,): i > 1) is a 
Poisson point process on (0, oa)x U with intensity ^“^d^v(di 7 ). 

The functions {/^; s e A) are known as spectral functions of 77 and (fT 2 li as de Haan’s spectral representation. When 
V is a probability measure, one can view Y{s) = /j as random variables on the probability space (U, U, v) and then (fT 2 li 
becomes (| 4 |i. Conversely, any representation (fT^ can be cast in the form (| 4 |i with a change of variables. Depending on 
the context one representation may be more convenient than the other. In terms of (fT2l i. the concurrence probability 
formula in (fTOl) becomes 


p{S\_,...,Sk) 



fsfu) 

max —— 
M,...,kf,{u) 


-1 


v(dM) 


v(di7), 


(13) 


and the proof is essentially the same. 
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One could expect from Definition |6] and Theorem |2] that the extremal concurrence probability depends on the 
distribution of the spectral process T in (|4]i or the choice of spectral functions in (fTSl i. These representations are not 
unique, but we will see in the theorem below that the extremal concurrence probability p{si,.. ■,Sk) depends only on 
the distribution of the max-stable process rj and not on the choice of the specific spectral representation. 

Theorem 3. For s\,..., Sk & X., k>2, we have 

k 

p{su...,Sk)^^{-iy ^ Ei^[logP;,(77(i^)<?7(ij),;6y}], (14) 

r=\ yc(l,...,it) 

\J\=r 


where f] is an independent copy ofrj. In particular when k — 2, 

p{si,S2) ^2 + Efj [logP;,{?7(ij) < fi(sj)J =1,2}]. 
Proof. Starting from (fTTI) and applying the inclusion-exclusion formula, we have 

p{su- ..,Sk)^^f, [a ({/(sj) > disj) for some ; 6 {1,..., k}])] 

k 


= ^fj 


2 2 (-l)-'A({/(.,) 6 7}) 

\J\=r 


(15) 


Since the cumulative distribution function of rj is P;, (?7(sj) < fjisj), j ^ J I = exp [-A ({/(s^) > i){sj), j 6 yj)], we get 

k 

p{su...,Sk)^'^(-iy ^ Ei^[logP;,(77(sj)<?7(sj),7 6 7 |]. 

/■=! yc|i,...,q 

\J\=r 


The simplification when A: = 2 is straightforward because when |7| = 1 we have 

Ej^ [logP;, {?7(ij) < iiisj),] 6 yj] = Ej^ {-'7('Syr*| = -1> 

as fjisjY'^ is a standard exponential random variable. □ 

In the remaining part of this section, we investigate some properties of the extremal concurrence probabilities. 
Surprisingly, although the two notions are different, we encounter strong similarities with the extremal coefficient ®. 
We recall that the extremal coefficient 9{s\, sf) takes values in [1,2], the lower and upper bounds correspond to perfect 
dependence and independence respectively. The next proposition states a similar result for the extremal concurrence 
probability. 

Proposition 1. For all si, sy 6 X, we have 

i) p(si, sf) — 0 if and only ifTj(si) and rjisf) are independent; 

ii) p(si, sf) — 1 if and only ifTj(si) and rjisf) are almost surely equal. 

The proof uses the following generalization and improvement of the upper bound in (|9]). 

Lemma 1. For all s\,..., Sk & X, k >2, we have p{s\, ..., s^) < E {minj=i_..j: T(sj)|. 

Proof. In the context of Theorem|3 we have (by conditioning on Y) 

Eyi max ^ ) > max y(s,)“*Ey fT(s/)| = | min T(s;) 

since Ey{?(iy)) = 1. This, in view of (ITOl i implies the desired result. □ 
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Proof of PropositionU] Equation (|9l) implies that p{si,S 2 ) = 0 if and only if 6 {si,S 2 ) = 2 which is equivalent to 
the independence of rjisi) and 77 ( 52 )- When 72 ( 51 , 52 ) = 1, Lemma [T] entails F(5i) = ^( 52 ) almost surely so that 
^(■si) = almost surely. It is easy to prove the converse implication: if 77 ( 51 ) and 77 ( 52 ) are almost surely equal, the 
same holds for E( 5 i) and E( 52 ) so that 72 ( 51 , 52 ) = 1 . □ 

Interestingly 72 ( 51 ,..., 5 <.) can be expressed via the extremal coefficients of another max-stable process. 

Proposition 2. Let p be a simple max-stable process as defined in (l4]i and ^,^ 1 , 772 ,... independent copies of it. 
Consider the simple max-stable process 


^( 5 ) = rnaxf;—-, 5 6^Y, 

i>\ pi(s) 

then 

k 

72(51,..., 5^) = ^(-1)'+' Yj 

|•=l yc|i,...,/t| 

\J\=r 

where 9^(sj, j € J) — — logP{^( 5 y) < 1,7 6 7). In particular 72 ( 51 , 52 ) = 2 — 07 ^( 51 , 52 ). 

Proof Clearly ^ is a simple max-stable process since both Y and 77 are non negative and E{F( 5 )/? 7 ( 5 )) = 1 for all s e X. 
We have 


Ejj [logP^ (77(57) < 77(57), 7 6 7}] 



max 


fliSj) I 


logPf {^(57) <1,767} = -^(57,7 6 7). 


The next corollary lists some properties of the extremal concurrence probability fun ction that closely par a llel th ose 
of the extremal coe fhcient function. In view of Proposition the proof follows as in ISchlather and TawnI 020031] or 
ICoolev et al.l II 2 OO 6 I 1 . 

Corollary. Let p: h p(o, h) be an extremal concurrence probability function associated to a stationary max-stable 
process in X for some arbitrary origin o € X and h 6 X. Then the following assertions hold. 

i) The function hi-^ p{h) is positive semidefinite; 

ii) The function h 1 —> p{h) is not differentiable at the origin unless p(h) — 1 for all h 6 X; 

Hi) If d > 1 and ifp is isotropic, then h 1 —» p(h) has at most a jump at the origin and is continuous elsewhere; 

iv) {2- p(hi H- h 2 )} < {2 - 72(/7i)){2 - p{h 2 )]forall hi,h 2 6 X; 

v) [2 - p{hi + h 2 )f < {2 - 72 (/ 7 i))" H- {2 - p{h 2 )f - 1 for all hi,h 2 eXandO<a<l; 

vi) [2 — p{h\ + /72))“ > {2 — p{h\)Y H- {2 — p(/i 2 ))“ ~ 1 for all hi, 112 6 X and a < 0. 

We conclude this section with an unexpected result that relates the bivariate extremal concurrence probability with 
the well known Kendall’s r. 

Theorem 4. For any max-stable process 77 , we have 72 ( 51 , 52 ) = t where r = E [ 57 ^ 77 ( 77 ( 51 ) — 77 ,( 5 i)) 57 ^ 77 { 77 ( 52 ) — 77 ,( 52 ))] 
is the Kendall’s t 0 /( 77 ( 51 ), 77 ( 52 )) and 77 , is an independent copy ofp. 

Proof. Let W - F{p{si), 77 ( 52 )) where F is the bivariate cumulative distribution function of ( 77 ( 51 ), 77 ( 52 )). Lrom ([TST i 
we have 72 ( 51 , 52 ) = 2 -H E(log W). But since ( 77 ( 51 ), 77 ( 52 )) is a bivari ate max-stable random vector, we know that 

P(W < w) - w-(l - t)w log w, 0 < w < 1 IIGhoudi et al.l 1199811 and hence, after some simple calculations, 

72(51,52) =T. □ 
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3 Formulas for extremal concurrence probabilities 


In this section we gather formulas for the extremal concurrence probabilities for some popular models of max-stable 
random vectors and processes. As we will see, it is not always possible to get explicit formulas, and in such situations, 
we propose to use Monte-Carlo methods. 


3.1 Closed forms 


Example 1 (Logistic model). The concurrence probability for the A:-variate logistic model, i.e., with cumulative dis¬ 
tribution 


F(zu---,Zk) = exp 


Z 



0 < or < 1, Zi,..., Zi > 0, 


is p(ii, ...,Sk)^ n;=i'(i - O'/;)- 


Recall that for this model independence is reached when a - \ while perfect dependence occurs as a i 0 and, as 
expected, for such situations we have p{s\,.. .,Sk) -Q and p{s\,... ,Sk) - 1 respectively. 


Proof. It is not difficult to see that the multivariate logistic model corresponds to the case where T in dUi is a pure 
noise process with margins such that P{y(s) < y) = exp[-{r(l - a)yY^^°’'\ where T is the Gamma function. Using 
Theorem|2] we have 


p{Si,...,Sk) =Ey 


Ey < max 


nsi) 


j=i,...,k Yisj) 


•U 


= Ey([-l0gU{?(5i),...,y(5,))f‘) 




1/Q'I 


T(k - a) 
T(k)T(l-a) 


1=1 


where the fourth equality used the fact that Yj%i Yisj) is a Gamma random variable with scale r(l - and 
shape k for which negative moments are known. □ 


Example 2 (Max-linear model). Consider the max-linear model pis) - max^^i.where Zi,... ,Z„ are 
independent unit Frechet random variables and some functions {i £ Cq, = !,...,«) such that - 

1 for all s e X. We have 


i) The concurrence probability equals 


p(su---,sk) = 'y^peisu---,Sk), 

e=\ 


pe{su---,sk) = 




. k (ff(sj) 


\m=\ 


with the convention that 0/0 = 0, a/0 = oo if a>0, and l/oo = 0. 


ii) The probability that component { dominates at sites si,.. .,Skis given by the term pf in ( fTbl l. i.e., 

pc{si,...,sk) = P{77(s^) = <pe{sj)Z(, ;■ = 1,.. .,k} 


(16) 


(17) 


Proof. Part i) is an immediate consequence of ( fT3l l. Indeed, let 1/ = {1,...,«) be equipped with the counting measure 
v{l) = ■ ■ ■ = v{n) = 1. By taking fs{m) - if mis), m 6 U, we obtain that the max-linear model has the representa¬ 
tion (fT2l i. The integral expression of the concurrence probability in (fOl) then becomes a sum of the terms p( in (fTbl l. 

Part ii) shows an intriguing fact that the concurrence probability p{si,.. .,Sk) for the max-linear model is the sum 
of the probabilities that one of the ii components dominates the rest. Indeed, by the max-stability property and the 
independence of the unit Frechet random variables Z,„’s, we have that the right-hand side of ( fTTb equals 

f ‘finis i) 1 

P < max max - —Z,„ <Z(>- P(flZi < Z 2 ), 

[ mi^e j=\....,k (fitiSj) ) 

where a = ‘fmisfiftisf). Equation (fTTl) follows from the fact that P(flZi < Z 2 ) = (1 - 1 - «)“*, a > 0 . □ 
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Example 3 (Chentsov random fields). Suppose that the p rocess Ff^) = where A c is a random set. Then, 
by analogy with the theory of symmetric a-stable process OSamorodnitskv and Taaaulll994l Chap. 8], the max-stable 
process 

? 7 (s) = max^,TA.(s), seX, 

i>i ' 

where T, s 1^, are independent copies of T s 1^, will be referred to as a Chentsov max-stable random field on X. 

For a Chentsov-type max-stable process we have 


p{Sx,...,Sk) = P({ii,.. C A I {ii, . . n A 0), 


(18) 


or less formally that the extremal concurrence probability is the conditional probability that the sites s\,.. .,Sk are 
covered by the random set A given that at least one of the sites is covered. 

Proof. The result is an immediate consequence of Theorem|2] Since the outer expectation in (fTOl i can be restricted to 
the event {miny=i. . ^ Y{sj) > 0), we have 


p(Su...,Sk) =Ey 

= Ey 


Er 

1a(u)\ 

[j=i,....k 1a(s;) J 

Ey. 



- -1 \ 
l|minj=i.t U(Sj)>0) 

-1 \ 
min lA(i;) 


E(min^=i,...,i: IaC^;)} 
E (max;=i,...,i IaC^))} 


(19) 


where in the second relation we used the fact that l{minj,i j, utip) - niin 7 =i,...,it lA('Sj) and that the event min,1 a('S/) > 
0 is equivalent to {lA(ii) = 1, ■. ■, lAi^k) =1)- This proves (flSl l. □ 

Example 4 (Extremal processes). The max-stable process {pis): s e [0,1]) is an extremal process if it has stationary 
and independent max-increments, i.e.. 


ipisi),p(sk)} = [iiZi, max{iiZi, (i 2 - si)Z 2 },..., max{iiZi,. ..,(sk- Sk-i)Zk}], 


where 0 < si < ■ ■ ■ < Sk and Z\,.. .,Zk are independent unit Frechet random variables. It can be shown that 

//(i) = max^,l[£/,A](i), i6[0,1], 

1>1 


where t/,’s are independent t/(0,1) random variables. Using our previous result on Chentsov random fields, we have 
for all 0 < Si < ■ ■ ■ < < 1 

^ P({si,...,s^}c [U,l]) ^ P(U< si) ^ s^ 

pKsu...,Sk) ^ 0 ) V{U<Sk) Sk’ 


where U ~ 11(0, 1). This result is not surprising since for this simple case, extremes are concurrent at locations 
0 < Si < ■ ■ ■ < Sk < 1 if p{s) has no jumps in the interval [si, s^,]. Hence using the independence and stationarity 
of the max-increments, the probability of the latter event is PjsiZi > (sk - si)Z 2 ) = si/sk, where Zi and Z 2 are two 
independent unit Frechet variables. 

Example 5 (Indicator moving maxima). In the context of (fT2li . if fsiu) = 1a,(m), for some sequence of measurable 
deterministic sets As, by using (fTTt . we obtain as in ( fT9l l that 


piSl, ...,Sk) 


v(ny=i,...,;tAj.) 

v(Uj=i,...aA^,)' 


( 20 ) 


















In the simple case /s(m) = 1a(m - ■?), i-e. = i + A with some deterministic set A, where v is the Lebesgue measure 
on W^, (I 20 I 1 implies 


p{s, s + li) — p(h) — 


|A n (h + A)| 
|A U {h + A)| 


CA(h) 

2\A\-CA(hy 


where |A| denotes the ^^-dimensional volume of A and caW - |A n (h + A)|. The latter function and hence the extremal 
concurrence probability function p{h) can then be obtained in closed form for many different sets. For example, in the 
case 77 is isotropic, i.e., A = {i 6 : ||jc|| < r) is the centered ball of radius r > 0 in Euclidean space, using the formula 

for the volume of the cap, we obtain 


CA(l|h||) = Cdr‘'B{d+i)/2,i/2 


||/ 2 ||( 2 r-||/ 7 ||)] 
2 r 2 j ’ 


Cd^ 




ril+d/iy 


where Ba^h(x) - B(a, b) ^ 11 “ *(1 “ 'd« is the cumulative distribution function of a Beta(a, b) random variable. 


3.2 Monte-Carlo methods 


It may happen that for some parametric max-stable models, explicit forms for extremal concurrence probabilities are 
not available but hopefully it is often possible to use Monte-Carlo methods to approximate the theoretical extremal 
concurrence probabilities with arbitrary precision. A naive strategy would consist in using (fTOt to devise a Monte- 
Carlo estimator, but it is wiser to take advantage of the closed forms of max-stable processes cumulative distributions, 
i.e., 

v{T](sj) < Zj,; = 1,.. .,k} = exp{-Fi,,...,,j(zi,. ..,z;i)), zi,.. .,z^ > 0, 

where is an homogeneous function of order -1. Rewriting (fTOl) . we found 

p{si,...,Sk) = Ey^[-log ( 77 ( 57 ) < Y{sj),j = l,...,k}] 

-Ey({y,,....,jy(5i),...,y(5,)]r') ( 2 i) 

which can easily be estimated by sampling independent copies of Y and computing the sample mean. We can often 
make use of antithetic variables to get more precise estimates. Note that specihc choice of the spectral process Y can 
lead to better strategies as we will illustrate in the following examples. 


Exampl e 6 (Brown-Resnick mo del). Let 77 be a Brown-Resnick stationary random held on X driven by a Gaussian 
process I Kabluchko et al. . 2009l] . That is, the processes y, in (HJi are equal in distribution to 


y(s) = exp{iy(s) - 7 ( 5 )), s e X, 


( 22 ) 


where W is a zero mean Gaussian random held with stationary increments and semi-variogram 7 , i.e. 2y{h) - 
e{w(/ 7)2) = E[{iy(i -I- h) - iy(5))2], s,h 6 A. 

For this model, the bivariate extremal concurrence probability function is given by 

77 ( 0 ,/ 7 ) = e([ 0(Z) +exp { 7 (/ 7 )- V2^z}(I){V27W-z}]^‘), (23) 


where Z ~ N(0, 1) has the standard normal distribution with cumulative distribution function <1). As expected p(o) - 1 
and p{h) ^ 0 as ||/ 7 || ^ 00 provided that the semi-variogram is unbounded, i.e., j{h) ^ 00 as ||/ 7 || ^ 00 . 

Proof. Without loss of generality, for Brown-Resnick processes we can assume that in (|4]i we have Y{o) - 1 almost 
surely. The bivariate cumulative function is given by P{ 77 ( 0 ) < zi,rj{h) < Z 2 ) = exp{-y/,(zi,Z 2 )), Zi,Z 2 > 0 with 


V/i(Zl,Z2) = —<h 
Zl 


Vr(/*)/2 ■ 


1 , Z2 

,- log — 

y2j{h) Z\ 


+ —d) 

Z2 


Vr (/^)/2 ■ 


1 , Zl 

- - log — 

y2y{h) Z2 


Equation (ISTT i together with the fact that {y(o), y(/7)} has the same distribution as {l,e rWj entails 

77(o,/r)-E({y,[y(o), y(/r)]r') = E{(y4l, exp{V^Z- 7 (/ 7 )}!)■'}. 

Equation (|2^ follows after straightforward simplihcations. □ 
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A po pular special case of the the Brown-Resnick family of models is the moving maximum storm model intro¬ 
duced bv ISmith 1 1990l] known also as the Gaussian extremal process. Consider the spectral representation (fT^ . where 
U - and v is the Lebesgue measure. Taking fsiu) - ipi:{s - u), where (f-^ is the multivariate Normal density with 
zero mean and covariance matrix S, we obtain the max-stable process 


77 ( 5 ) = max - Ui), s € X. 

i>l 


(24) 


Then, the following are true; 

i) The process 77 belongs to the family of degenerate Brown-Resnick models in (l22l i with lT(i) = s eW^ and 

ii) Consequently, the concurrence probability function of 77 is given by (l2Tt with y{h) - 

Proof. Let 77 be a Brown-Resnick process with variogram 2y{h) - that is, we can assume without loss of 

generality that in the spectral characterization we have Y{s) = exp(i^Z“'Z - with Z ~ A(0, Z). Then for all 

,..., i/t e X and zi,. ■ ■,Zk > 0 we have 


-logP{ 77 (s;) < ZjJ = 1,. 


. .,k] = Ez I max z• * exp(sJZ *Z-sJZ 


= (27r) ^^^|Z| r max z-'exp(-^y’^Z 'y-HiJZ 

= I max zyVsCs;-y)dy- 

jRdj=l,...,k 


The last relation equals the negative log cumulative distribution function of the moving maxima in (l24l i. □ 

Example 7 (Schlather and extremal-f processes). Let 77 be an extremal-f process on X, i.e., the processes T, in (01 are 
equal in distribution to 


Y(x) = Cy max{0, IT(i))'', Cy 



seX, 


where v > 1 and VT is a stationary standard Gaussian process with correlation function p. The Schlather process is 
obtained when v = 1 . 

The corresponding extremal concurrence probability function is 


77(0, h ) = 


Ty^iiT) + {p(h) + crihyTrTyyi 


pjh) 

o-(h) 


o-(h)(p(h) + o-(h)T) j 


\ 

1 {p{h)+cr(h)T >0) 


where o-(h) - - p(hf}l(l + v) and T is a Student random variable with v - 1 - 1 degrees of freedom and cumulative 

distribution function Ty+i. 

Proof. For the notational convenience, we write shortly p - p{h), cr - cr{h). To obtain the desired result, we use a 
different spectral representation. Comparing the two cumulative distribution functions, one can show that 

{77(o),77(/7)) = max^,{y,(o), Yiih)] 


with Yi, i > 1, i.i.d. copies of the bivariate random vector 


{YiofYih)} 


|2(1, max{0,p - 1 - crTY), with probability 1/2 
| 2 ( 0 , c), with probability 1/2 


10 














with c = 1 - E[max(0,p + crTY] such that E(F) = (1,1). Equation (|2TI) yields 


p{o,h) = E([y,{?(o), = E([y,.{l,(p + o-TY)]-^ l(p+<.r>0|) , 


with ODavison et al.Ll2012 l 


Z\ 



1 


1 1T J 

1 P 

1 



-H — 

— 

f H- Ty+l 1 


+ — 

— i 

1 ^ 

cr 

Ui/ j 

[ 22 1 

1 ^ 

cr 

U 2 / j 


and straightforward simplifications give the announced result. 


□ 


4 Statistical inference and asymptotic properties 


4.1 Sample concurrence probability estimators 

In this section we define a sample concurrence probability estimator by blocking the data and study its basic properties 
as well as the optimal choice of the block-size. We conclude with a methodological improvement of the estimator 
based on permutation bootstrap. 

Let Xi = {Xi{sj): j - I,... ,k}, i = 1,... ,n, he random vectors in k>2. Partition the data into non-overlapping 
blocks of size m < n, and define the sample concurrence probability estimator 


Pm — Pmi.^1 5 ■ - • 5 


j [n/m] 

[n/m] z 


(25) 


The max statistic above is simply an indicator of whether or not we have concurrence in the r-th block, i.e., whether 
one of the vectors dominates the componentwise maximum of the rest in the r-th block of size m. 

Assuming that X\,.. .,X„ are independent and identically distributed, the above estimator is the sample mean of 
[n/m] independent Bernoulli(p,„) random variables, where /?„, is as in © with n replaced by m. Therefore, 

in f \ Pm (1 Pin ) 

E(p„,) = Pm, Var(/7,„) = ————, 

[n/mj 

i.e., pm is a unbiased estimator for p„,. 

As argued in the introduction, a major drawback of the sample concurrence probability /?,„ is that it depends on the 
sample size m and it is thus more sensible to focus on the limiting extremal concurrence probability p - p{s \,..., sY). 
The sample concurrence probability estimator is biased for p with mean squared error 

MSE(p„,) = {pm -pf+ . (26) 

[n/mJ 

Although the bias term Pm-p is difficult to estimate in general, for the max-stable case a precise expression is available. 

Proposition 3. Assume that X is a max-stable process, then pm is non-increasing in m and satisfies 

A ^ ^ (1 -P) 

0 < Pm - P < - , m > 1 . 

m 

Furthermore, if p < 1, then there exists an integer re {I,.. .,k — 1) and a positive constant Cr 6 (0,1 — p\ such that 
iPm - p) ~ Crim'' asm ^ oo. 

The proof expresses {pm-p) via the distribution of the extremal hitting scenario and is postponed to Apr)endix lA.2l 
The above result suggests that an asymptotically optimal choice of the block size can be made to minimize the rate 
of the mean squared error in (l26l l. In view of Proposition[3] for the general case 0 < p < 1, 

MSE(pm) = f—f + oim-^'') + 0(m/n) 

\ m' / n 
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Taking the derivative with respect to m we see that the optimal rate corresponds to 2rc^m ' ~ p(l - p)n as 
n —> oo. Hence the block size that asymptotically minimizes the mean squared error is 


'«MSE(n) ~ 


Irc^n 

P(l-P) 


l/(2r+l) 


n —> oo. 


(27) 


For this rate-optimal mean squared error we obtain MSE(p„,) oc n 2 ''/( 2 r+i) 

Remark. The constants r, Cr and p in ( l27b are unknown. The precise expressions for r and Cr involve multiple 
concurrence probabilities, i.e., when two or more events contribute to the maximum. In principle, pilot estimates of 
these parameters could be obtained and used as plug-ins in dZTl i. Since the cases r > 2 correspond to very specific 
dependence structures, in practice, we recommend using the conservative choice r = 1 and = 1. 

The following result establishes the asymptotic behavior of the sample concurrence probability estimator. The 
proof is given in Appendix lA.SI 

Theorem 5. Suppose that 0 < p < 1 and let m — m{n) be such that njinin) —> oo and —> d e (0, oo] as 

n — > oo. Then 

s/n/mipm-p )—> -p)|, n^oo 

where l/oo is interpreted as zero. 


Remark. In Theorem |3 we encounter a typical tradeoff between rate-optimality and bias. In particular, the mean 
squared error optimal choice of m as in ( IZTT i corresponding to d = {2rCrf'l^^''^^\p{\ - yields the limit 

distribution 


N 


Vf(i -p) 


2r 


p{\-p) 


which has non-zero mean. On the other hand, the rate sub-optimal choices where A - oo yield Normal unbiased limits. 

In the case p = 0, we have no optimal block size or optimal rate estimation. Still, the following asymptotic result 
can be useful. 


Theorem 6. Suppose that p = 0 and let m — m{n) be such that nlm{n) oo and —> d e (0, oo] as 

n ^ oo. If A < oo, then 

(nlm)p,„ —> Poisson(cr/d''^'), n —> oo. 

Otherwise, if A = oo, then P(p,„ = 0) —» 1 as u —» oo and hence a„p,„ —» 0 in probability as n ^ oo for any sequence 

an > 0 . 


Remark. When k = 2 it is possible to get an unbiased estimator for p based on a slight modification of p„,. Indeed for 
this specific case, implies that p = {mpm — l)/(m - 1) and hence the estimator 


Pm ~ 


mp,n - 1 

OT — 1 


(28) 


is unbiased for p as p„, is unbiased for p„,. 

We conclude this subsection with a brief methodological improvement of the sample concurrence probability 
estimator p,„ based on permutation bootstrap. The idea is to compute the estimator p,„ for several independent random 
permutations of the sample Xi,... ,X„. Then the average of the resulting estimator would have a lower variance and 
the same mean p^. 

Formally, this procedure is justified by the following simple observation based on Rao-Blackwellization. Consider 
the lexicographic linear order in R*’, denoted <, and let A(i) < A( 2 ) < ■ ■ ■ < A(„) be the sorted sample obtained from 
X\,.. .,Xn. The independence of the A,’s and the continuity of their marginals entails that the above ordering is strict 
with probability one. 

Let T{X\,... ,Xn] - (2f(i),..., A(„)). It can be shown that T is a sufficient statistic for the parameter p,„ = 
Pniisi,.. .,Sk) and the Rao-Blackwell theorem implies the following propostion. 
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Proposition 4. For p*^ = E(p,„ | T) we have E (p*„) = /?„, and 


^[ip:n-P,nf}<E[{p,„-p,„f}. 


Moreover we have 


Pm ~ I Pm{^iT{\)^ ■ ■ ■ ^Xcr(n)], 
n\ 


crG5„ 


and where S„ denotes the set of all permutations of{\, 
An alternative expression for p*^ is 



where di = 2Li ^iXk<Xii and = 0 ifdf < m - 1. 


(29) 

(30) 


(31) 


Proof Let/(x; p,„), x 6 be the density of X = {2f(ii),... ,X(i(i;)) with respect to some dominating measure A. The 
fact that for the likelihood, we have 


L{p,„-,Xi,i = 1,,..,«) = ]~[/(Xi;p„) = j~[/(X(,■);;>,„), 

1=1 1=1 

shows that T - {^^(i),... is a sufficient statistic for p,„. The inequality in ( |29] | follows by appealing to the 
Rao-Blackwell Theorem or simply applying the conditional form of Jensen’s inequality. 

The independence of the Xi’s and the lack of ties (with probability one) in the lexicographic order imply 

P{Xa-(i) = X(,), / = 1 ,..., n} = 

n\ 

for all cr e S„. This shows that p* - E{p,„ \ T) is expressed as in (l30l l. 

By the definition of the sample concurrence probability estimator, we get 


1 1 [n/ffi] 

-y — y 1 

nl ^ [n/m] 


{sample concuiTence occurs within (Xo.(,+(r-i);n) : / = 1,... ,m)| 


creS, 


(■) ^ 


ample concuiTence occurs within (Xi : i G S)} 


S£P,n(n) 


where Pinin) is the collection of all subsets S c {1,... ,n) of m elements. Given a subset S 6 Vm{n) and io e S, it is 
easy to see that dominates (Jf,),es if and only if 5 \ {I'o) is included in the set {k - 1,... ,n',Xic < We deduce 
that for a given index /q, the number of subset S e Vm{n) such that dominates (Jf,),e 5 is equal to J- Formula ( |3T] ) 
follows easily. □ 

The above result shows that the estimator p*^ is superior to in terms of mean squared error. From a numerical 
point of view, formula (ISTT i is much more computationally efficient than (l30t . We shall refer to p*, as to the sample 
concurrence probability bootstrap estimator. It is significantly better, in practice, than the simple sample concurrence 
probability estimator p„, and therefore in applications we recommend using only p*„. As indicated above, in the case 
of bivariate concurrence, the bias of p,„ can be removed. As in Relation (l28l l. we obtain the following unbiased 
modification of for the case of pairwise concurrence 


P 


m 


rnPm - 1 

m — 1 


The respective performances of p,„, p*„ and p*„ are analyzed in Section|5] 


(32) 
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4.2 Extremal concurrence probability estimators 

Although the sample concurrence probability can be easily estimated, deriving an estimator for the extremal concur¬ 
rence probability seems at first sight difficult since in (fTOt the stochastic processes Y and Y are not observable. Recall 
that for statistical purposes we often assume that the pointwise block maxima, e.g. pointwise annual maxima, are 
distributed according to some max-stable process and thus we observe 77 but not Y. Fortunately, Theorem[3]enables us 
to estimate p{si,.. .,Sk) without the need of observing Y. 

Based on rp,.. independent copies of 77 , one possible estimator for pisi,.. .,Sk) is to consider the sample 
counterpart of (fT4l) . i.e.. 


p(Si,---,Sk) 


k ^ n \ M n 

Z<-»' Z 

r=\ (=1 k=\ 

\A=r 


\nt(Sj)<r]i(Sj), jeJ] 


(33) 


In the innermost summation, we include the index k - i to ensure that the logarithm is always well defined. 
Although estimator dS?! seems natural it has undesirable properties since it is not linear in the data and is thus likely to 
show some bias for small sample sizes. Also the study of its asymptotic properties seems delicate. However, to reduce 
the bias of this estimator, it is always possible to use a Jackknife procedure. 

Fortunately, in the bivariate case, it is possible to get an unbiased estimator based on Theorem 01 This theorem 
suggests the simple estimator p = p(si, S 2 ) - t, where f is the Kendall’s t statistic, i.e.. 


p = T ■ 


77(77 — 1 ) 


^ sign{77,(si) - 77/si)}sign{77,(s2) - 77 /S 2 )}, 


(34) 


l<i<j<n 


that is well known to be unbiased and such that I Dengler , 2010[ Theorem 4.3] 

V^(t - r) ^ N(0, ctI), crl - 15 Var[T{ 77 (si), 77(52)) - ^{77(51)) - ^{77(52))]. 


Although the asymptotic variance cr^ is hard to evaluate as it requir es knowledge of the dependence structure, in 
practice it can be accurately and consistently estimated using Jacknife llSchemt)eilll987l] . 


4.3 Integrated concurrence probabilities and area of concurrence cell 


As we will see in Section | 6 ] the above methodology can be used to provide bivariate concurrence probability maps 
s p{so, s) centered at a given location sq. Such maps show how fast the dependence in extremes decreases when 
moving away from 5o- A drawback of this approach is that one may produce one such map for every choice of an 
origin so and the choice of an origin is hence quite arbitrary. To bypass this issue, we propose to consider the integrated 
concurrence probability 


7 (so) = 



so 6 


Intuitively, this quantity measures how fast the dependence in extremes decreases when moving away from sq . In¬ 
terestingly, it can be related to the notion of concurrence cell and its area. Consider the Poisson process representation 
of the max-stable random field 77 = ( 77 ( 5 ): s e A) in Q. Recall that we have a concurrence of extremes at sites so and 
s if 77 (so) = 0(so) and 77 ( 5 ) = 0(s), for the same (p e Let C(so) denotes the random set of all sites s that are in a 
concurrence relation with so. This set will be referred to as the concurrence cell containing the site sq. 


Proposition 5. For any site so e X, we have /(sq) = E{|C(so)|) where |C(so)| is the d-dimensional volume of C{so)- 

Proof Observe that the concurrence probability satisfies p{so, s) = E { 1 c(jo)('*) 1 ^^d that the volume of the concurrence 
cell is given by 

|C(so)| = r lc(i„)(i)ds. 

Jx 

The result follows by applying the Tonelli-Fubini’s theorem. □ 
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Block size Block size 

Figure 1: Evolution of the root mean squared error for p„, (left) and pl, (right) as the block size m and the sample size 
n increase. These estimates were obtained from 2000 Monte-Carlo samples sampled from a Brown-Resnick model 
with semivariogram y{li) — hf 1.621. This semivariogram was chosen such that the theoretical extremal concurrence 
probability is p - 0.5. The red circles indicate the optimal block sizes as defined by dZTl i and their corresponding 
optimal root mean squared error (|26]) . 


We will provide and discuss in Section |6] some maps of the integrated concurrence probability so liso) that 
allow to evaluate at each location so the dependence in extremes around so. For a detailed study of the properties of 
the concurrence cells associated to a max-stable rando m field and of the tessellation of t he entire domain generated by 
the concurrence cells, please refer to the recent work of iDombrv and Kabluchkol 020141] . 


5 Simulation study 


In this section, we analyze the performance of the pairwise sample concurrence probability estimators p,„, p*, and p*, 
defined in (l25l l. (OTT i and (l32l i respectively, and that of their extremal counterpart p in (l34l i. Since the latter estimator 
relies on the max-stability assumption while the former three assume that observations belong to the max-domain of 
attraction, we need to handle both situations. The fi rst one is well known and consists in sampling from max-stable 
processes using the methodology of ISchlathen 0200211 . In the second situation, to be able to control the degree to which 
the model differs from a max-stable one, we consider the following partial maxima 


ijis) = — max *T,(i), s e X, 

no !=1,...,K0 


(35) 


where T, are as in (IHi, U\,..., Una independent t/(0, 1) random variables and for some suitable no 6 N. By construc¬ 
tion, fj belongs to the max-domain of attraction of in (HI and in so me sense can be viewed as a truncation of the 
spectral representation in (|4]i (see, e.g., the proof of Proposition 3.1 in IStoev and TaqquI 1120051] .') 

We first focus only on the sample concurrence probability estimators, i.e., p„, and p*,, and analyze their perfor¬ 
mance with respect to the block size m and the sample size n. Based on a Monte-Carlo simulation, Figure[T]shows the 
evolution of the root mean squared error as the block size grows. As expected, both estimators become increasingly 
more efficient as the sample size grows and, as seen from ( l29b , the permutation estimator p*, is more efficient than 
p ,„—independently of the block size m and the sample size n. The circles on the plot indicate the asymptotically 
optimal block size in (l27l i. which are valid olnly for max-stable data. As expected the observed optimal block sizes 
are in good agreement with the theoretical ones. In practice, however, since the data are not exactly max-stable, we 
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Figure 2: Evolution of the root mean squared error for p as the theoretical extremal concurrence probability p and the 
number of spectral function no in dlSl l increase. These estimates were obtained from 2000 Monte-Carlo samples of 
size n with, from left to right, n - 25,50,100,500. 


recommend using slightly larger values of m so as to ensure that the block-maxima are closer to a max-stable model 
but also to take into account that data usually exhibit serial dependence, e.g., daily observations. 

We now investigate the performance of the extremal concurrence estimator p. Figure |2] shows the evolution of 
the root mean squared error as the number of spectral functions no in (l35] l and the theoretical extremal concurrence 
probability increase. As expected, as the sample size n grows, the estimator p becomes much more efficient. Interest¬ 
ingly, for small sample sizes, p appears to be fairly robust to the lack of max-stability in the data, i.e., no < This 
is not true anymore for larger sample sizes since, as expected, p becomes increasingly more efficient as the number of 
spectral functions increases. 

Finally, we compare the performance of the sample concurrence probability estimators p*^ and with their 
extremal concurrence counterpart p. To compare the two types of estimators on a fair basis, we analyze their behaviour 
when the simulated data are either perfectly max-stable or non max-stable, but in the domain of attraction of a max- 
stable distribution. 

Figure|3]shows boxplots of the sample p*,, unbiased sample p*, and extremal concurrence probability estimators p, 
based on 2000 Monte-Carlo realizations of both a Brown-Resnick and extremal-f models. Recall that we focus here on 
pairwise concurrence probabilities. In this case, the extremal concurrence probability coincides with Kendall’s t and 
therefore, the estimator p in (l5^ is in fact unbiased for the case of max-stable data. This is confirmed by the results in 
Figure 13 As expected, the variability of all estimators decreases as the sample size grows; the extremal concurrence 
probability estimator being the most precise one. Since the simulated data are max-stable, we can see that the sample 
concurrence probability estimator is biased even when the sample size is large while the remaining two estimators are, 
as expected, unbiased. Overall the extremal concurrence probability appears to be the best estimator provided that the 
data are max-stable. 

To corroborate this finding, Table [T]reports Monte-Carlo sample means and standard deviations of these estimators 
as the assumption of max-stability becomes more accurate, i.e., as the number no of spectral functions in (ITST i grows. 
As expected, when the max-stability assumption is most unreasonable, i.e., no = 1, all estimators show a substantial 
bias with the extremal concurrence probability estimator p having the largest bias while the unbiased sample one 
p*, the lowest. As the assumption of max-stability becomes increasingly more accurate, the bias of the unbiased 
sample concurrence and extremal concurrence estimators improve. When this assumption holds exactly (indicated 
by no = oo), both of these estimators exhibit essentially no bias as stipulated by the theory and seen in Figure |3] 
The sample concurrence probability estimator appears to be biased in all situations—the bias being less significant 
as the number of spectral functions is larger. Interestingly, whatever the estimator considered, the bias and variance 
appear to increase as the theoretical extremal concurrence probability value p becomes smaller. Overall the extremal 
concurrence probability estimator p in (l34l i has the lowest variability. 
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Figure 3: Boxplots of the sample (red/ middle), unbiased sample (green /right) and extremal (grey/left) concurrence 
probability estimators at distance lags h — 1,2,3,4. The boxplots were obtained from 2000 independent estimates. 
From left to right: the sample size is respectively 25, 50, 100 and 500. The top panel corresponds to an extremal-t 
model with v = 5, and correlation function p(li) = exp(-/i/10). The bottom panel corresponds to a Brown-Resnick 
model with semi variogram y{h) - hl3. For each panel, the solid line represents the corresponding theoretical extremal 
concurrence probability function. 


6 Concurrence of temperature extremes in continental USA 


In this section, we apply the developed methodology to estimate the probabilities of concurrence associated with 
extreme temperatures—both extreme cold and hot events. The data consists of daily temperature minima and maxima 
recorded at 424 weather stations over the period 1911-2010. The spatial distribution of the se station s is giv en in 
Figured This data set, as a subset of the United States Historical Climatological Network lUSHCNI ll2014ll . was 
chosen as it meets very high data quality standards and involves fewer than 2.4% missing values while spanning the 
entire territory of continental US. It can be freely downloaded from http://cdiac.ornl.gov/ 

To avoid any seasonal influence on our results we decided to analyze minima and maxima for each season sepa- 



o 



Figure 4: Left: Spatial distribution of the 424 weather stations. The triangles indicate the selected stations for the 
analysis — upward: daily maxima, downward: daily minima. Right: The seasonal extrema time series of the selected 
stations. 
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Table 1: Performance of the sample (pl,), unbiased sample (p*„) and extremal (p) concurrence probability estimators. 
The table report the sample mean and the standard deviation in paren based on 2000 Monte-Carlo replicates. The 
data are either simulated from an extremal-t model with correlation function p(h) = exp(-/i/10) and v - 5 degrees of 
freedom or from its truncated representation with no extremal functions. Throughout this simulation study the block 
size is held fixed tom — 10, independently of the sample size n. 




p = 0.25 



p = 0.50 



p = 0.75 



Pm 

Pi 

P 

Pm 

Pi, 

P 

Pi, 

Pi, 

P 





Sample 

size n = 20 





no = 1 

0.41 (0.24) 

0.35 (0.26) 

0.47 (0.13) 

0.64 (0.22) 

0.60 (0.24) 

0.71 (0.09) 

0.83 (0.15) 

0.81 (0.17) 

0.87 (0.05) 

no = 10 

0.34 (0.24) 

0.27 (0.25) 

0.31 (0.14) 

0.57 (0.23) 

0.52 (0.26) 

0.58(0.12) 

0.79 (0.17) 

0.77 (0.19) 

0.80 (0.07) 

no = 15 

0.33 (0.24) 

0.27 (0.25) 

0.30 (0.15) 

0.56 (0.23) 

0.52 (0.26) 

0.56 (0.12) 

0.78 (0.17) 

0.76(0.19) 

0.78 (0.07) 

no = 00 

0.33 (0.24) 

0.27 (0.25) 

0.25 (0.15) 

0.55 (0.24) 

0.50 (0.26) 

0.50 (0.13) 

0.77 (0.18) 

0.75 (0.20) 

0.75 (0.08) 





Sample 

size n = 50 





no = 1 

0.41 (0.13) 

0.35 (0.14) 

0.47 (0.08) 

0.65 (0.10) 

0.61 (0.12) 

0.71 (0.05) 

0.84 (0.07) 

0.82 (0.07) 

0.87 (0.03) 

no = 10 

0.34 (0.13) 

0.26(0.14) 

0.31 (0.09) 

0.57 (0.12) 

0.52 (0.13) 

0.57 (0.07) 

0.79 (0.08) 

0.76 (0.09) 

0.80 (0.04) 

no = 15 

0.33(0.13) 

0.25 (0.14) 

0.29 (0.09) 

0.56 (0.12) 

0.51 (0.13) 

0.56 (0.07) 

0.78 (0.08) 

0.76 (0.09) 

0.79 (0.04) 

no = CO 

0.32 (0.13) 

0.25 (0.14) 

0.24 (0.09) 

0.54 (0.12) 

0.49 (0.14) 

0.50 (0.08) 

0.77 (0.09) 

0.74 (0.09) 

0.75 (0.05) 





Sample ; 

size 11 = 100 





no = 1 

0.41 (0.08) 

0.35 (0.09) 

0.46 (0.06) 

0.65 (0.07) 

0.61 (0.07) 

0.71 (0.03) 

0.83 (0.04) 

0.82 (0.04) 

0.87 (0.02) 

no = 10 

0.34 (0.09) 

0.26 (0.10) 

0.31 (0.06) 

0.57 (0.08) 

0.52 (0.09) 

0.57 (0.05) 

0.78 (0.05) 

0.76 (0.05) 

0.80 (0.03) 

no = 15 

0.33 (0.09) 

0.26(0.10) 

0.29 (0.06) 

0.56 (0.08) 

0.51 (0.09) 

0.55 (0.05) 

0.78 (0.05) 

0.76 (0.06) 

0.78 (0.03) 

no = CO 

0.33 (0.09) 

0.25 (0.10) 

0.25 (0.07) 

0.55 (0.08) 

0.50 (0.09) 

0.50 (0.05) 

0.78 (0.05) 

0.75 (0.06) 

0.75 (0.03) 


rately. We focus on the concurrence of extreme cold (minima) during the Fall and Winter seasons—generally color- 
coded in blue; and extreme hot (maxima) during the Spring and Summer seasons—generally color-coded in red. The 
right panel of Figure |4] shows the times series of these seasonal extrema for four selected weather stations. These 
stations were selected as they recorded the top seasonal records over the whole spatial and temporal domains. We can 
see that all four time series of seasonal extremes (cold in blue and hot in red) at these stations appear to be stationary 
without any clear tem poral trend. This is i n contrast with the generally accepted trend of about 0.2°C per decade for 


average temperatures OStocker et al.L 1201311 . 


Figure |5]plots the estimated spatial distribution of the extremal concurrence probabilities function for each season, 
relative to the chosen station. More precisely, for a given origin location sq, the maps display estimates of the pairwise 
concurrence probability p(sq, s) as a function of s. These maps were obtained by first computing the estimator (|34]i 
over all 423 pairs of stations (^n, 5 ) a nd then interpolated using thin plate splines (on logit scale) provided by the R 
package fields I Nvchka et al.L 2014|] . As expected, the highest concurrence probability occurs in the neighbourhood 
of the selected stations independently of the season. The areal extent of high concurrence probabilities, however, seem 
to be larger for minimum temperatures (cold extremes) than for maximum temperatures (hot extremes). This finding 
is consistent with the physical notion of entropy, i.e., when the ambient temperature is higher (Spring and Summer 
seasons), the entropy is greater and hence involves less spatial dependence than for cooler temperatures leading to 
smaller probability of simultaneous extremes. This difference can be also attributed to the fact that extreme cold 
temperatures are often due to high-pressure systems, which tend to linger longer and cover a larger spatial area than 
warm fronts giving rise to concurrence of extreme hot events. 

Although Figure|5]displays interesting patterns, it has the drawback of being dependent on the choice of the origin, 
i.e., the selected station. As stated in Section |3 it is possible to bypass this hurdle by considering the area of concur¬ 
rence cell. Figure |6] plots the estimated spatial distribution of the concurrence cell area for the preindustrial period, 
i.e., 1910-1975, and the postindustrial one, i.e., 1976-2010. To emphasize the possible impact of anthropogenic in¬ 
fluences, the ratio of these two cell areas is also reported. We can see that during the last sixty years the expected cell 
area for winter minima have increased of about 30% over the whole USA while there is a decrease of about the same 
amount for summer maxima. These findings indicates that today’s climate shows cold spells that have a larger impact 
than in th e beginning of the 20th century while hot spells are more localized. Our results agree with the conclusions 
drawn by Field et al.l 020121] who states that “there is evidence from observations gathered since 1950 of change in 
some extremes”. These changes in the concurrence patterns of summer extremes can be attributed to global warm¬ 
ing since an increase in entropy generally leads to more “mixing” in the system and hence less dependence leading 
to smaller areas of concurrence. The changes in concurrence patterns of extreme cold events, however, are harder 
to explain. They may be triggered by structural changes in important climatological mechanisms such as the Arctic 
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Figure 5: Maps of the extremal concurrence probability for the four selected stations. Top left: Fall (September, 
October November), top right: Winter (December, January, February), bottom left: Spring (March, April, May) and 
bottom right: Summer (June, July, August). 



Figure 6: Estimated spatial distribution of the expected extremal concurrence cell areas—in squared degree, i.e., 
around 1000 km~. From left to right: 1910-1950, 1951-2010, and their ratio (1951-2010 at the numerator). Top: 
Winter minima, bottom: Summer maxima. 
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Figure 7: Spatial distribution of the estimated concurrence cell areas anomalies in squared degree, i.e., around 
1000 km^ for winter minima (top) and summer maxima (bottom). The data were stratified into three classes: El 
Nino, La Nina and the base class “La Nada”. The left panels show the anomalies for El Nino, the right ones La Nina. 


Oscillation. 

Finally, we consider another cut of the data by stratifying according to an important climate phenomenon known 
as the El Nino Southern Oscillation (ENSO). Positive ENSO (El Nino) refers to the event of a warm-up of the surface 
water in the central and east-central equatorial Pacihc ocean. It is well known that years with high ENSO have a 
general warming effect in North America during the winter season. The opposite effect of negative ENSO (La Nina) 
is characterized by a cool-down in the same area of the Pacihc and it generally l eads to unusua lly cold winters in 
the northwestern part of the US, northern California and the north-central states GrahamI 1 1999l] . Eigure [T] shows 
estimates of concurrence cell areas anomalies for winter minima and summer maxima. These anomalies were dehned 
as pointwise deviations from the expected cell area obtained from La Nada seasons, i.e., neither El Nino nor La Nina 
seasons. We can see that La Nina does not seem to have an impact on the spatial coverage of winter minima but that 
El Nino seems to induce more massive cold extremes over the whole USA. Eor the summer season. La Nina seems to 
reduce the spatial extent of heat waves over the whole USA while El Nino has a less pronounced impact—although it 
generally yields larger spatial coverages especially along the East coast. 


7 Discussion 

In this paper we introduced a new framework for the analysis of dependence of extremes: the extremal/sample concur¬ 
rence probability. This tool plays a similar role to that of the extremal coefficient but has the benefit, as a probability, 
of being more interpretable and intuitive. Theoretical properties and closed forms of these concurrence probabilities 
have been established and several estimators have been proposed. A simulation study has shown that the proposed 
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estimators work well in practice and that they give a new insight about the dependence of extremes, such as the spatial 
distribution of the expected concurrence cell area of extreme temperature in the continental US. 
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A Proofs 

A.l Proof of the continuity of 11 in Theorem [1] 

Consider a sequence T'„ ^ T* in ^p{[Q, oo]^ \ (0)) and let T = {i/r,-: i > 1) and = {i^„j : i > 1). As compact subsets 
of [0, oo]* \ (0) are bounded away from 0, we can choose e > 0 such that 


mnxi>n//i{sj) > E, 

if/iisj) E, j - . ,k, i > \. 


Then Kg - [0, oo]*’ \ [0, e)^ is a compact set, T' n /fg has finitely many points i/^i,..., and no point of T' lies on 
the boundary dK^. The convergence T'„ ^ T' entails that for n large enough, T'„ n Kf. has the same number of points 
,..., ilJn,N^ that can be reordered in such a way that i/r„^ t/r, as « —» oo, f = 1,..., Ns- 

Assume that the maxima maxi<,<Ar^ j - 1,..., k, are uniquely attained. Then the hitting scenario n(T') 

is well defined and depends only on {i/r, : / = 1,... ,Ns]- By the convergence —» i/r, for large n the maxima 

maxi<,<A^^ ifr„j(sj) are uniquely attained so that the hitting scenario n('T„) is well defined and depends only on {tf/nj: i = 
1,... Ae). It is not difficult to see (although tedious to write formally) that the convergence : i - 1,..., As) —> 
(tk,: i - 1,..., Ns] implies n('P„) = n('I0 for large n. This proves the announced continuity for the mapping 11. 

A.2 Proof of Proposition |3] 

We shall prove below the following formula, which may be of independent interest. 

Lemma 2. In the context of Proposition]^ we have 


k 



(36) 


where n is the extremal hitting scenario and \n\ its number of components. 

Proposition [3] follows directly from (|36]l. Indeed, the monotonicity of p,„ is immediate and since p = P(|7r| = 1), 
we have 



If p < 1, then at least one of the probabilities P(|7r| — I), £ - 2,... ,kis non-zero and by ( |3^ the asymptotic equivalence 
iPm — p) ~ Crim'', OT —> oo holds where r e [I,.. .,k- l)is the smallest integer, such that Cr = P(|7r| = r H- 1) > 0. 

Proof of Lemma^ Since Zi,... ,Z„, are independent with the same distribution as rj, we can suppose from dD that 


Zi(s) - maxf(s), s e X, 
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with Oi,..., <[)„, independent copies of O. By max-stability, tj - m ' maxi<,<„, Z, has the same distribution as 77 and 
<i) = (ft 6 d)/) has the same distribution as d). We consider the following events 

A - {sample concurrence occurs for Z\,.. Z,„) 

B - {extremal concurrence occurs for 77 ) 

Clearly, P(A) = and P(A) = P(A n B) H- P(A n B'’) with the complementary set of B. We analyze the two terms 
separately. 

Observe first that if extremal concurrence occurs then we also have sample concurrence. Indeed, if B occurs, 
then one function of d) dominates all the others at (ii,..., Sk). This function is of the form with (f) 6 d),, for 
some 7 = l,...,77i, showing that Z, dominates Zi,...,Z,„, i.e., we have sample concurrence. Hence B c A and 
P(A n B) = P(B) = p. We now consider the second term and the event A n B'^, i.e., sample concurrence occurs in 
Zi,..., Z,„ but not extremal concurrence for 77 . Let n be the hitting scenario of 77 . We know that B‘^ is equivalent to 
\n\ > 2, i.e. the maximum at locations (si,..., Sk) is attained by at least two functions in d). These functions are of the 
form m~^(f)j, 1 < j < (, with i - | 7 r| and (pj e d), for some 1 < ij < m. If A is also realized, i.e., some Z, dominates 
Zi,..., Z,„, then we must have i\ - ■ ■ ■ - ic - i. Note, however, that since the point processes d>,, 7 = 1,..., ??i are 
independent and identically distributed, any given function e d> s U'” j 7 « 'd>,, independently from the others, 
has equal chance of coming from any one of the m point processes 7 « 'd>,, i - I,.. .,m. Therefore, the probability that 
all { functions contributing to the maximum at sites (ii, ..., s^) are assigned to component 7 is Since there are m 
possible choices for the index 7 , we deduce 


k 

P(A n B") = ^ P(|7f| = 
{=2 


Equation < |3^ follows. □ 

A.3 Proofs of Theorems |5] and |6] 

In the context of these two theorems, we have 


Pm - p ~ cpn \ in 00 , (37) 

for some r e {l,...,k-l) and > 0—cf. Proposition [3] Let = [nlm]pm - where are iid 

Bernoulli(77,„). 

Proof of Theorem\^ Let B„ = Var(5'„) = {nlrn\p,„{\ — p,„) and introduce the cumulative distribution function 


F„{x) = P [b;'/2{5^ _ E{S„)] < = pi 


[?7/77i] 


I Pni (f Pm) 

The Berry-Essen theorem (see e.g. Theorem V.2.3 in Petrov 1 1975ll l implies that 


(Pm - Pm) < ■* L xeR. 


sup \F„(x) - d)(x)| < AL,„ 

xeR 


(38) 


where A is an absolute constant, d>(x) denotes the standard Normal cumulative distribution function and 


[n/m] 

i=l 

Using that E(|^/„ - Pm\^) ^ Pm(^ - Pm) and straightforward algebra, we obtain 

L„ < [{n/in)p^{l - Pm)T'^^ ■ (39) 


22 










Since 0 < /? < 1, we have /?m(l - /?„,) ~ p(l - p) > 0 as « —> oo and, for an arbitrary choice of m - m{n) - o(n) as 
n — » oo, we have L„ —> 0. This, in view of (l38l l yields 


■ Pm) -> N{0,p(l- p)}. 


(40) 


By using (iTTl i and Slutsky’s theorem, we obtain the final result. □ 

Proof of Theorem^ In case d < oo, since p -0, from (iJTl i. we have p,„ ~ Crim'', n oa, and hence 

[n/m]p„, ~ n ^ oo. (41) 


Thus, the standard Poisson convergence for the Binomial(/7m, [n/m]) random variables S„ — [n/m]pi„ yields the result. 
When d = oo, by (HTt . we have x„ = [nlm\pm — > 0 as n — » oo. Therefore, 

/ \ [«/'"] 

P(A« = 0 ) = (1 - ^ ~ In/^) ~ ®^P(“-^«) —^ 1’ u ^ oo, 

which completes the proof. □ 
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